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Preface 



These are notes for part of a course on advanced quantum mechanics given to 4th 
year physics students. The only prerequisites, however, are a basic knowledge of the 
Schrodinger and Heisenberg pictures of standard quantum mechanics (as well as the will- 
ingness to occasionally and momentarily suspend disbelief). Thus the material could 
easily be presented at an earlier stage. I covered the material in five 3- "hour" lec- 
tures (1 "hour" = 45 minutes) and this time constraint (there are other topics that I 
wanted to cover as well in the course) dictated the level of detail (or lack thereof) of 
the presentation. 

One of the aims of these lectures was to set the stage for a future course on quantum field 
theory. To a certain extent this motivated the choice of topics covered in these notes 
(e.g. generating functional are discussed, while concrete applications of path integrals 
to non-trivial quantum mechanics problems are not). 

These notes do not include an introductory section on motivations, history, etc. - such 
things are best done orally anyway. My own point of view is that the path intgral 
approach to quantum theories is simultaneously more intuitive, more fundamental, and 
more flexible than the standard operator - state description, but I do not intend to get 
into an argument about this. Objectively, the strongest points in favour of the path 
integral appoach are that 

• unlike the usual Hamiltonian approach the path integral approach provides a man- 
ifestly Lorentz covariant quantisation of classical Lorentz invariant field theories; 

• it is perfectly adapted to perturbative expansions and the derivation of the Feyn- 
man rules of a quantum field theory; 

• it allows in a rather straightforward manner to calculate certain non-perturbative 
contributions (i.e. effects that cannot be seen to any order in perturbation theory) 
to S-matrix elements. 

The motivation for writing these notes was that I found the typical treatment of quantum 
mechanics path integrals in a quantum field theory text to be too brief to be digestible 
(there are some exceptions), while monographs on path integrals are usually far too 
detailed to allow one to get anywhere in a reasonable amount of time. 

I have not provided any referenes to either the original or secondary literature since 
most of the material covered in these notes is completely standard and can be found in 
many places (the exception perhaps being the Gelfand-Yaglom formula for fluctuation 
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determinants for which some references to the secondary literature are given in section 
3.5). 

No attempt at mathematical rigour (not even the pretense of an attempt) is made in 
these notes. 

Updated corrected or expanded versions of these notes will be available at 

http : //www.blau. itp .unibe . ch/Lecturenotes .html 

If you find any mistakes, or if you have any other comments on these notes, complaints, 
(constructive) criticism, or also if you just happen to find them useful, please let me 
know. 
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1 The Evolution Kernel 



1.1 Review: The time-evolution Operator 

The dynamical information about quantum mechanics is contained in the matrix ele- 
ments of the time-evolution operator U(tf,U). For a time-independent Hamiltonian H 
one has 

U(t f ,t i ) = e-( i / h H t f- t *)H (L1) 

whereas the general expression for a time-dependent Hamiltonian involves the time- 
ordered exponential 

t%, w -r ■ d.2) 

It satisfies the evolution equation 

ih£-U(t f ,U) = H(tf)U(tf,t i ) , (1.3) 

with initial condition 

U(t f = ti,U)=l . (1.4) 

As a consequence, 

nt f ) = u{t f , umu) (i.5) 

satisfies the time-dependent Schrodinger equation 

ih df f ^f) = H(tfMtf) (1-6) 

with initial condition 

*(ti) = V'(ii) • (1-7) 
Another key-property of the time-evolution operator is 

U(t f , U) = U(t f ,t)U(t, U) U<t<t f . (1.8) 

In words: evolution from time U to tf is the same as evolving from t{ to t, followed by 
evolution from t to tf. 

1.2 Definition of the Propagator / Kernel 

Our main interest will be in matrix elements of U in the position representation. We 
will denote these by (the letter K stands for "kernel") 

K(x f ,Xi;tf-ti)=<x f \U(t f ,ti)\xi> . (1.9) 
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This can also be interpreted as the transition amplitude 

K(xf,Xf,tf - ti) =< Xf,tf\xi,U > , (1-10) 
where (in the time-independent case) 

\ Xhti>=e m)uH\ Xi> (L11) 

etc. Note that this is not the Schrodinger time-evolution of the state \xi> - this would 
have the opposite sign in the exponent (this is related to the fact that U(tf,ti) satisfies 
the Schrodinger equation with respect to tf, not U). Rather, this state is characterised 
by the fact that it is the eigenstate of the Heisenberg picture operator xn(t) at t = ti, 

x H (t)\x,t>=e( i / h ) t ^xe-^/ h ) tA e^/ n ) tA \x>=x\x,t> . (1.12) 

Another way of saying this, or introducing the states \x,t >, is the following: In general 
the operators xu{t) and £#(£') do not commute for t / t' . Hence they cannot be 
simultaneously diagonalised. For any given t, however, one can choose a basis in which 
xn(t) is diagonal. This is the basis {\x,t >}. 

Matrix elements of the evolution operator between other states, say the transition am- 
plitude between an initial state \ipi > and a final state \ipf >, are determined by the 
kernel through 

<ipf\U(tf,ti)\ipi> = J dxf J dxi < ipf\xf X Xf\U(tf,U)\xi >< Xi\ipi > 

= J dxf j dxi i(>f(xf)ipi(xi)K(xf,Xi;tf - U) . (1-13) 

1.3 Basic Properties of the Kernel 

Here are some of the basic properties of the Kernel: 

1. With properly normalised position eigenstates, one has 

lim < Xf, tf\xi, U >=< Xf\xi >= S(xf — Xi) . (1-14) 

2. The kernel satisfies the Schrodinger equation with respect to (tf,Xf), i.e. 

[ihd tf - H(xf,pf = (h/i)d Xf ,tf)] < Xf,tf\xi,U >= . (1-15) 
More generally, for any ket \<p > the wave function 

^^{x,t)=<x,t\4>> (1.16) 
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is a solution of the Schrodinger equation, with 

*</,(x, t) = J dxi < x, t\xi, U X Xi,ti\<f) > 

= J dx i K(x,x i ,t-t i ) i If l p(xi,ti) . (1-17) 

3. Since we can restrict our attention to evolution forwards in time, one frequently 
also considers the causal propagator or retarded propagator 

K r (xf,Xi;tf - U) = Q(tf - ti)K(xf,Xi;t f - U) (1-18) 

where 6 is the Heavyside step function, @(x) = 1 for x > and @(x) = for 
x < 0. It follows from the above two properties of the kernel, and O'(x) = 8(x), 
that the retarded propagator satisfies 

[ifrd tf -H(xf,pf = (h/i)d Xf ,tf)]K r (xf,Xi;tf-ti) = ihS(tf-U)S(xf-Xi) . (1.19) 

Thus the retarded propagator is a Green's function for the Schrodinger equation. 

4. A key property of the kernel is the convolution property 

/+0O 
dx K(xf,x;tf — t)K(x, xf,t — U) (1-20) 
-00 

(for t arbitrary subject to the condition tf > t > tj). This follows from the 
property (1.8) of the time-evolution operator by inserting a complete state of 
states in the form 

/oo 
dx \x >< x\ . (1.21) 
-oo 

One then finds 

< Xf,tf\xi,U > = <Xf\U(tf,t)U(t,U)\xi> 

/oo 
dx < Xf\U(tf,t)\x >< x\U(t,U)\xi > 
-oo 

/oo 
dx < x f,tf\x,t >< x,t\xi,t{ > , (1-22) 
-oo 

which is the claimed result. 

5. Finally, we can also expand the kernel in terms of energy eigenstates 

ip n (x) =< x\n > (1.23) 

as 

K(xf,Xi;tf-ti) = <x f \e-( i / n Wf- ti ^\ Xi> 

= ]T <x / |e-(^)(*/-^|n><n|x i > 

n 

= Y,*~^ ,mf ~ U)En MxfW n {xi) ■ (1-24) 

n 
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1.4 Example: The Free Particle 



For the free particle, with Hamiltonian H = p 2 /2m, it is straightforward to determine 
the kernel K = Kq and verify explicitly all the above properties. By inserting a complete 
set of momentum eigenstates \p >, 

< x \ p >= ^ =e m) P x (L25) 

V27r?i 

(i.e. by Fourier transform), one finds 

K (x f , Xl ;t f -U) = <x / |e-(^)(*/-^ 2 /2m| Xi> 

= f + °° dp < x f \p > e "WW*/ " ti)P 2 / 2m < p \ Xl > 

J — oo 

= — [ + °° dvettWte&f ~ Xi ^ ~ (*/ -*i)P 2 / 2 H (126) 

27T/ii_ 00 ^ • v • 7 

Using the Fresnel integral formulae from Appendix A, one thus finds 

^''^''-v-J^mr^" (L27) 

Note that the exponent has the interpretation as the "classical action" , i.e. as the action 
So of the free particle evaluated on the classical path x c (t) satisfying the free equations 
of motion and the boundary conditions x c (tf^) = Xfj, 

S [x c ] = m {X /- Xl) * . (1.28) 
L J 2 (tf-U) 

This at this stage rather mysterious fact has a very natural explanation from the path 
integral point of view. 

1.5 Exercises 

1. Verify the results (1.27) and (1.28) for the evolution kernel of the free particle. 

2. Verify that (1.27) is normalised as in (1.14), and satisfies the convolution property 
(1.20) and the free particle Schrodinger equation (1.15). 

3. Prove the basic integral identity (A. 31), 

/ + oo °°^e-^V2 = ( det _i)- 1/2 . (L29) 

for Gaussian integrals (A a b is a real symmetric positive matrix) , and calculate the 
integral 

fV*e-^^V2 + Ja* a (1 . 30) 

J — oo 

by completing the square (remember that A a ^ is invertible). 
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2 Towards the Path Integral Representation of the Kernel 



2.1 From the Kernel to the Short-Time Kernel and Back 

In general, it is a difficult (if not impossible) task to find a closed form expression for the 
kernel. However, the convolution property allows us to reduce the determination of the 
finite-time kernel to that of the short-time (or even infinitesimal time) kernel K(x,y;e) 
and, as we will see later on, this allows us to make some progress. 

First of all we note that we can write (once again in the time-independent case, but the 
argument works in general) 

< z f \e-WWts ~ U)H lx . >=< Xf] (^-m^Hy^ > (21) 

We think of this as dividing the time- interval [U,tf] into N equal time-intervals [t k , ifc+i] 
of length e, 

e = t k+1 -t k = (t f -U)/N . (2.2) 

Here k = 0, . . . , N — 1 and we identify tf = tjv and U = to- We can now insert N — 1 
resolutions of unity at times t k , k = 1, . . . ,N — 1 into the above expression for the kernel 
to find 

N-l N-l 

K(xf,Xi;tf - U) = [ Yl / dx k ][ Yl K(x k+1 ,x k ;e = t k+1 - t k )} , (2.3) 
k=i k=0 

where Xf = xn and Xi = xq. 

Now this expression holds for any N. But for finite N, the kernel is still difficult to 
calculate. As we will see below, things simplify in the limit e — > 0, equivalently iV — > oo. 
In this limit, the kernel 

N-l N-l 

K(x f ,x i ;t f -h)= \im [H dx k }[ JJ K(x k+1 , x k ; e = (t f - U)/N)} (2.4) 

k =i k=0 

is determined by the short-time kernel K(x k+ i,x k ; e) as e — > 0. 

2.2 The Trotter Product Formula and the Dirac Short-Time Kernel 

Formally, it is not difficult to see that in the limit N — > oo we only need to know 
this short-time kernel K(x k+ ±, x k ; e) to order e. If everything in sight were commuting 
(instead of operators) , the argument for this would be the following: one writes 

e x = (e x / N Y = (1 + x/N + 0(1/N 2 )) N (2.5) 
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and compares with the formula 



e x = lim (l + x/N) N (2.6) 

N— >OD 

to conclude that in the limit N — > oo the subleading 0(N~ 2 ) terms in (2.5) can indeed 
be dropped. 1 Of course, to establish an analogous result for (unbounded) operators 
requires some functional analysis. 

We will now assume that the Hamiltonian is of the standard form 

2 

H = T(p) + V(x) = ^ + V(x) . (2.7) 

We are thus interested in determining the kernel 

K(x k+1 ,x k ;e)=<x k+1 \e-( i / h X f + V)\ Xk> . (2.8) 

Provided that we can justify dropping terms of 0(e 2 ), things simplify quite a bit. Indeed, 
even though T and V are non-commuting operators, eT and eV commute up to order 
e 2 , because their commutator is e 2 [T, V]. Thus, using the Baker-Campbell-Hausdorff 
formula one has 

e -(i/f l )e(f + V) =e -(i/h)ef e -(i/h)eV + 0{e 2 ) (2g) 

To justify dropping these 0(e 2 ) commutator terms, however, one needs some control 
over the operator [T, V] which should not become too singular. 

Concretely, what one needs is the validity of the Trotter product formula 

e A + B = ^A/N + B/Ny li^ e A/N e B/ N y {2AQ) 

for A = f and B = V. 

This identity is not difficult to prove for bounded operators. The case of interest, 
unbounded operators, is trickier. The identity holds, for instance, on the common 
domain of A and B, provided that both are self-adjoint operators that are bounded 
from below. Once again, we will gloss over these functional analysis complications and 
proceed with the assumption that it is legitimate to drop the commutator terms (while 
keeping in mind that this assumption is not valid e.g. for the Coulomb potential!). 



lr The identity (2.6) can be proved directly, e.g. via a binomial expansion or by showing that 

4- lim (1 + x/N) N = lim (1 + x/N) N 

aX N—>ao JV->oo 
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With the above assumptions, to order e we can write the short-time kernel as 

K(x k+U x k ;e) =< x k+1 \e e -(.WeV^ > . (2 . n) 

To diagonalise the operator T we introduce a complete set of momentum eigenstates 
\p k > with 

< x \ Pk >= ^ =e m) Pk x (212) 

and calculate 

K(x fc+1 ,x fc ;e) = f + d Pk <x k+1 \e-^/ h ^\p k X Pk \e-^/ n ^\x k > 

J — CO 

2 

= f + °°dp k <x fc+1 |p fc ><p fc |x Jfc >e-( i /»Mfe + ^ fc )) 

J— oo 

= — r dp.e^^^f^ - H (P^ (2.13) 

This is a lovely result, first obtained by Dirac in 1933. In the exponential, where we once 
had operators, we now encounter the Legendre transform of the classical Hamiltonian, 
i.e. the Lagrangian! Indeed, if we identify 

Xfc+i - x k _ Xfc+i - x k _^ dxk_ ^ 
e tk+i — tk dt 

as a discretised time-derivative, the exponent takes the classical form p k x k — H(p k ,x k ). 
We can implement the Legendre transformation explicitly by performing the Gaussian 
(Fresnel) integral over p k (see Appendix A) to find 

K(x k+1 ,x k ;e) = -L [ + °° dp k e (^)M^ + i - **) - + V(x k ))] 

m e {i/h)eL(x k ,x k ) (215) 
2irihe v ' ; 

where the Lagrangian L is 

L(x k , x k ) = m(x ™- Xk? - V{x k ) ^ m f- V (x k ) . (2.16) 

2.3 The Path Integral Representation of the Kernel 

Having obtained the above explicit formula for the short-time kernel in terms of the 
Lagrangian, we can no go back to (2.4) to obtain an expression for the finite-time 
kernel. We can use either the phase space expression (2.13) or the configuration space 
expression (2.15). This iteration of Dirac's result is due to Feynman (1942). 
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In this way we arrive at 

N-l m N-l 



*(»„,;.,-« = lim (H f *"Alti f ^Si^ - *<»,*.)] 

Note that in the phase space representation of the kernel there is always one more 
momentum than position integral. This is a consequence of the fact that each short- 
time propagator contains one momentum integral whereas the position integrals are 
inserted between the short-time propagators and the two end-points X{ and xj are not 
integrated over. 

In the e — > or N — > oo limit, we interpret the points Xk as defining a continuous (but 
almost certainly nowhere differentiable - see below) curve, any two successive points 
being joined by a straight line. That is, we think of them as defining a curve x{t) with 
endpoints 

x(ti)= Xi x(t f ) = x f (2.18) 

and 

x k = x(U + ke) . (2.19) 
Likewise, we think of the p k as definig a curve p(t) in momentum space such that 

Po=p(U) p k =p(U + ke) . (2.20) 
With this interpretation the exponents in the integrand of the kernel can be written as 

N-l _ rtf 



lim e £ [ Pk Xk+1 Xk - H(p k , x k )\ = I 1 dt [p{t)±{t) - H(p(t), x(t))] 

N-l ft j 

lim e ]T L ( x k, x k ) = / dt L(x(t), x{t)) = S[x(t);t f , U] (2.21) 



In the same spirit, we formally now write the integrals as integrals over paths, intro- 
ducing the notation 

r x(t f )=x f 



/oo rx(tf)=Xf 

dx k ] = / D[x(t)} 
-oo Jx{U)=Xi 

j im i ff r = I D W)i*«*\ ■ ( 2 - 22 ) 



lim 

fe=l 
N-l 



With this notation, we can now write the kernel as a path integral 

r x(tf)=x f _ r _ f _ r _ _ f t 



Jx{ti)=Xi J 

M r {t f )=x f D[x ^ ]e (i/h)s[x(t)] _ (223) 

J x(ti}=Xi 
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This is an integral over all paths x{t) with the specified boundary conditions and, in 
the first version, all paths p(t) with t 6 [U,tf]. 

Some remarks: 

1. The first line, a phase space path integral, is valid for general Hamiltonians 
H(p,x) = T{p) + V(x) (possibly time-dependent). The measure appears to be 
an infinite-dimensional analogue of the canonical Liouville measure. The latter 
is of course invariant under canonical transformations. This should however not 
lead one to believe that the path integral representation of the kernel also enjoys 
this invariance. Indeed, this cannot possibly be true since it is well known that 
under a canonical transformation any Hamiltonian can be mapped to zero (the 
Hamilton-Jacobi transformation) and hence into any other Hamiltonian, while the 
kernel depends non-trivially on the Hamiltonian. 

2. To pass to the second line, a configuration space path integral, we used the explicit 
quadratic form T = p 2 /2m to perform the Gaussian integral over p. The derivation 
of the path integral for a Hamiltonian with a velocity dependent potential, such 
as the magnetic interaction (p — A) 2 /2m, is more subtle (the discretised version 
requires a "mid-point rule" for the Hamiltonian) and will not be discussed here. 

3. As mentioned above, as limits of piecewise linear continuous paths these paths are 
continuous but not differentiable. Indeed differentiability would require existence 
of a finite limit of (xjt+i — Xk)/t as e — > 0. But x^+i and x^ are independent 
variables, and hence there is no reason for the difference Xk+i — £fc to go to zero as 
e — > 0. Hence the paths entgering the above sum/integral are typically nowhere 
differentiable. Evidently, then, things like x(t) require some (perhaps stochastic 
or probabilistic) interpretation, but we will not open this Pandora's box. 

4. Finally, M is formally an infinite normalisation constant, 

/ mN \ N/2 

M= lim ^ (2.24) 

N->oo\2mh{tf -U) ) 1 ; 

(we have replaced e by (tj — U)/N) whose sole purpose in life is to make the 
combined expression J\f times the integral well defined, finite and equal to the left 
hand side (provided that all our functional analysis assumptions are satisfied). 

We could have absorbed much of the prefactor into the definition of the measure by 
writing (2.17) as 

- « - (^) 1/2 ^jn £ (^r-^ws 

(2.25) 
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and denning 

rx(t f )=Xf _ N ~} r°° / m X 1 / 2 

/ D\x{t)}= lim [TT / — dx k ] . (2.26) 

Then the kernel is 

Another useful normalisation of the measure is, as we will see in section 3, such that 
M [ X{tf) = Xf DixWeWWW)] = I jn f*f)=°f 6[x(t)]e (i/R)S[x(t)] . 

Mu)=xi V 2 ^Htf - U) Mu)=xi 

(2.28) 

We will obtain a more informative expresion for J\f, and hence the relation between D[x] 
and D[x], involving the determinant of a differential operator, in the next section. 



2.4 The Path Integral Representation of the Partition Function 
In analogy with the quantum statistical (or thermal) partition function 

Z(P) = Tre~P H (2.29) 

one defines the quantum mechanical partition function, or the partition function for 
short, as the trace of the time evolution operator, 

Z(t f ,U) = Tr U(t f ,U) . (2.30) 

In the time-independent case, this becomes 

Z(t f ,U) = Z(t f -U) =Tre - (Vft)(*f ~U)H ? (2 31) 

and thus the quantum mechanical and thermal partition function are formally related 
by continuation of the time interval (tf — U) to the imaginary value 

t f -ti = -ih/3 . (2.32) 

Evaluating the trace in a basis of energy eingenstates |n >, one finds 

Z{t f -t i ) = Y J ^~^ i/K ^ t} ~ ti)En ■ (2-33) 

n 

On the other hand, evaluting the trace in a basis of position eigenstates \x >, one 
obtains 

/+oo f+OO 
dx < x\U(tf,U)\x >= / dx K(x,x;t f -U) . (2.34) 
-OO J — CO 
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This means that in the discretised expression (2.17) for the kernel there are now an 
equal number of momentum and position integrals. In the continuum version (2.23), 
setting Xf = xi = x means that one is integrating not over all paths but over all closed 
paths (loops) at x, and the integration over x means that one is integrating over all 
closed loops, 



Since we now have an equal number of x- and p-integrals, in terms of the measure D[x] 
(2.26) the partition function reads 



Comparison with (2.33) shows that, if we are able to calculate this path integral over 
all closed loops we should (in the time-independent case only, of course) be able to read 
off the energy spectrum. 

The relation between statistical mechanics and quantum mechanics at imaginary time 
is rather deep. In particular, for quantum field theory this implies a relation between 
finite temperature quantum field theory in Minkowski space and quantum field theory in 
Euclidean space with one compact (Euclidean "time") direction. This has far-reaching 
consequences (none of which will be explored here). 

2.5 Back from the Path Integral to the Schrodinger Equation 

Since inside the path integral one is dealing with classical functions and functionals 
rather than with operators, fairly simple "classical" manipulations of the path integral 
can lead to non-trivial quantum mechanical identities. 

As an example, consider the "trivial" statement that the path integral is invariant under 
an overall shift 




(2.35) 




(2.36) 





(2.37) 



of the integration variable, 




which follows from the (presumed) translation invariance of the measure. 
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Infinitesimally, with y(t) = 5x(t) and 

S[x(t) + Sx(t)\ = S[x{t)} + SS[x(t)\ , (2.39) 
this statement reduces to 

r {tf)=Xf D[x{t)] 6S[x{t)] e (i/n)s[x(t)] = (240) 

J x(ti)=Xi 

Sometimes this statement is paraphrased as "the path integral of a total derivative is 
zero" , 

^'^(i)] « eWWLo . (2.41) 
i)=xi oa;(t) 

Consequences derived from such identities are known as "Schwinger-Dyson equations" 
in the quantum field theory context. 



rx(t- 
Jx{U) 



Since the variation of the action with Sx(t) vanishing at the endpoints gives the Euler- 
Lagrange equations, 

«W"i-jC'*(s-is) fa w- (2 ' 42) 

and since (2.40) holds for any such 5x(t), we deduce that the classical equation of motion 
are valid inside the path integral (!), 

cr^Cs-sS)-^ 1 - - <-» 

This is the path integral version of Ehrenfest's theorem. 

Since this was so easy, let us generalise this a bit and consider variations of the path 
which fix x(ti) but change x(tf), 

Sx(U) = Sx(t f ) / . (2.44) 

This means that we are calculating 

d 

5 < xf,tf\xi,U >= 5x(tf)-^— < xf,tf\xi,U > . (2-45) 
In this case, the variation of the action is 

z(t)+p(tf)5x{t f ) , (2.46) 





(dh 


d dL 


\dx 


dt dx 


r* 1 


( dL 


d dL 


\dx 


dt dx 
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where p(tf) = pf is the canonical momentum at t = tf. This leads to the standard 
Hamilton-Jacobi relation 

We deduce that 

A <x t | x ^> * ( X{tf)=Xf D[x{t)]p(t f )e^/ n ) S ^ , (2.48) 
oxf n J x (ti)=xi 

which is nothing other than the familiar statement that in the position representation 
one has 

In the same way, by considering a variation of tf, and using 

L{x f ,x f ) = ^S[x(t)] = ^ + Pf Xf (2.50) 

(modulo the Euler-Lagrange equations), and therefore 

^ = -H(x f ,p f ) , (2.51) 

one obtains a path integral of the statement (1.15) that the kernel satisfies the Schrodinger 
equation (Exercise), 

(ih^-H(x f ,p f = n ^)) [ X{tf)=Xf D[x(t)]em) S W)] =0 . (2 . 52) 

V % <>■'■/ J Jx(ti)=Xi 

We had initially derived the path integral from the standard description, and now 
we have employed path integral manipulations to pass back from the path integral to 
the Schrodinger representation. This essentially completes the (formal) proof of the 
equivalence of the path integral description of quantum mechanics and the standard 
Schrodinger representation. One piece of the dictionary that is still missing is how to 
translate matrix elements other than the basic transition amplitude < Xf,tf\xi,U > 
into the path integral language. This will be the subject of the next subsection. 



2.6 Convolution, Correlation Functions and Time-Ordered Products 

As we saw in section 1 and above, a crucial property of the kernel is the convolution 
property (1.20). How is this encoded in the path integral representation (2.23) 

K(x f , Xi - t f -U)=M [ X(tf)=Xf D[x(t)]e (*/ft)S[s(t); *fM (2.53) 

J x(ti)=Xi 
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of the kernel? We need to consider seperately the integrand and the measure. As far as 
the integrand is concerned, the action S[x(t); tf, U] obviously satisfies 



S[x(t);t f ,ti] =S[x(t);t f ,to] + S[x(t);t ,ti] (2.54) 
for any tf > to > U since 

rtr rtr rto 

/ (•••)= / (•••)+ / (•••) • (2-55) 

Jti Jto Jti 

And for the path integral measure we have 

rx(tf)=Xf f+oo rx(tf)=Xf rx(to)=xo 

/ D[x{t)} = dx D[x(t)] / D[x(t)} (2.56) 

Jx(ti)=Xi J — CO Jx(to)=XO J x(ti)=Xi 

In words: we can perform the integral over all paths from X{ to Xf, by considering 
all paths from Xi to xq and xq to xj for some fixed xo = x(to) and then integrating 
over xq. Taken together, the statements about the action and the measure imply 2 the 
convolution property (1.20). 

So far we have only discussed the transition amplitude < Xf,tf\xi,U >, but it is also 
possible to represent matrix elements of operators as path integrals. The most natu- 
ral operators to consider in the present context are products of Heisenberg operators 

XH(h)xH(t 2 ) ■ ■ 

We begin with a single operator xn(to) with tf > to > U and, to simplify the notation, 
we will from now on drop the subscript H on Heisenberg operators. By elementary 
manipulations we find 

/+oo 
dx(t ) < Xf,tf\x(to)\x(t ),t X x(t ),t \xi,ti > 
-oo 



/. 



'OO 

+oo 

dx(t ) < x f ,tf\x(t ),t >< x(t ),t \xi,ti > a; (t ) 



Turning this into a statement about path integrals, using the convolution property, we 
therefore conclude that 



<x / ,* / |x(t )|x i ,* i >= f X(tf) X/ Z)[x(t)]x(to)e(^)^W;i/,^ . ( 2 .57) 

Jx(ti)=Xi 



Thus the insertion of an operator in the usual prescription corresponds to the insertion 
of a classical function in the path integral formulation. While this is very charming, and 
in line with the replacement of the Hamiltonian operator by the Lagrange function and 



At least as long as one pretends that M = 1 or that Af has somehow been incorporated into the 
definition of a suitably regularised path integral. This is the recommended attitude at the present level 
of rigour (better: non-rigour), and one that we will adopt henceforth. 
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(2.58) 



its action, it also immediately raises a puzzle. Namely, since x(t±) and xfo) commute, 
does the insertion of x{t\)x{t2) into the path integral calculate the matrix elements of 
x{ti)x(t2) or x(t2)x(t\) or . . . ? To answer this question, we reverse the above calcu- 
lation, but this time with the insertion of x{t\)x(t2)- In order to use the convolution 
property of the kernel or path integral, we need to distinguish the two cases £2 > h and 
t2 <t±. Then one finds 

[ X(tf)=Xf D[x(t)] x(ti)a;(t 2 ) e {i/K)S[x{t)-t f ,U] 

J x(ti) = Xi 

= f X{tf)=Xf D [ x{t) ] X (t 2 ) x ( tl ) e m)S[x(t);t f , U) 

J x(ti) = Xi 

{ < Xf,tf\x(t2)x(ti)\Xi,ti > t 2 >h 
\ < Xf,tf\x(ti)x(t2)\Xi,ti> t 2 <ti 

Using the time-ordering operator we can summarise these results as 

[ X{tf) J Xf D[x(t)} x(h)x(t2) e (i/fr)S[x(ty,t f ,U] =< s/j t / |r( i r(ti)s(t2))|x i ,t i > . 

(2.59) 

This immediately generalises to 

[ X{tf) ~ Xf D[x(t)} x(h) . . . x(t n ) e (Vft)^(i); t f , U] =< Xfjtf \ T ^ tl) . . . x(t n ))\ Xi ,U > 

Jx(ti)=Xi 

(2.60) 

Thus the path integral always evaluates matrix elements of time-ordered products of 
operators. This explains why the ordering inside the path integral is irrelevant and how 
the ordering ambiguity is resolved on the operator side. 

As a final variation of this theme, we consider the path integral representation of tran- 
sition amplitudes between states other than the position eigenstates \x,t >. To obtain 
these, we simply use the formula (1.13), 

< ipf\U(tf,ti)\ipi >= JdxfJdxiifj* f (xf)ifji(xi)K(xf,Xi;tf-ti) , (2.61) 

and the path integral expression for the kernel K(xf,Xi;tf — U). 



2.7 Comments 



We have thus passed from a formulation of quantum mechanics based on the Hamilto- 
nian (and operators and Hilbert spaces) to a Lagrangian description in which there are 
only commuting objects, no operators. In this framework, the quantum nature, which 
in the usual Hamiltonian approach is reflected in the non-commutativity of operators, 
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arises because one is instructed to consider not only classical paths, i.e. extrema of the 
action (solutions to the classical equations of motion) but all possible paths, weighted 
by the exponential of (i/h) times the action. 

This provides an extremely intuitive picture of quantum mechanics in which one is 
never led to ask questions like "which path did the electron take?". One also sees very 
explicitly how classical mechanics emerges in the limit fi — > 0. In that case, the path 
integral will be dominated by contributions form the extrema of the action (by the usual 
stationary phase approximation), i.e. precisely by those paths that are solutions to the 
classical equations of motion (see the discussion in section 4.6). 

It is important to note that at this point the path integral notation that we have 
introduced is largely symbolic. It is a shorthand notation for the N — > oo limit (2.17). 
Provided that all our assumptions are satisfied, this is just another (albeit complicated 
looking) way of writing the propagator. 

However, the significance of introducing this symbolic notation should not be under- 
estimated. Indeed, if one always had to calculate path integrals as the limit of an 
infinite number of integrals, then path integrals might be conceptually interesting but 
that approach would hardly be an efficient calculational tool. One might like to draw an 
analogy here with Riemann integrals, defined as the limits of an infinite sum. In practice, 
of course, one does not calculate integrals that way. Rather, one can use that definition 
to establish certain basic properties of the resulting infinite sum, symbolically denoted 
by the continuum sum (integral) J, and then one determines definite and indefinite 
integrals directly, without resorting to the discretised description. Of course, in this 
process one may be glossing over several mathematical subtleties (which will ultimately 
lead to the development of measure theory, the Lebesgue integral etc.), but this does 
not mean that one cannot reliably calculate simple integrals without knowing about 
these things. 

The attitude regarding path integrals we will adopt in the following will be similar in 
spirit. We will deduce some properties of path integrals from their "discretised" version 
and then try to pass as quickly as possible to continuum integrals which will allow us to 
perform path integrals via one "functional integration" instead of an infinite number of 
ordinary integrations. Once again, this will be sweeping many important mathematical 
subtleties under the rug (not the least of which is "does something like what you have 
called D[x] exist at all?"), but that does not mean that we cannot trust the results that 
we have obtained. 3 

3 By the way, the answer to that questions is "no", but that is irrelevant - it is simply not the right 
question to ask. A more pertinent question might be "can one make sense of J D[x] exp(i/h)S[x] as 
something like a measure (or linear functional)?". 



21 



For some historical comments on 11 efforts to give an improved mathematical meaning to 
Feynman's path integral formulation of quantum mechanics" , see e.g. 

J.R. Klauder, The Feynman Path Integral: An Historical Slice, quant -ph/0303034. 
2.8 Exercises 

1. Using the Baker-Campbell-Hausdorff formula, determine the operator X, defined 



2. Generalise the derivation of the path integral to systems with d > 1 degrees of 
freedom. Assume that the Hamiltonian has the standard form 



where x = (x 1 , . . . , x d ) etc. 

3. Spell out the proof of (2.52) (the path integral satisfies the Schrodinger equation) 
in detail. 



by 



e -(i/h)e(f+V) =e -{t/h)ef e -{ l /h)eV e -(e/hfX 



(2.62) 



to order e. 




(2.63) 
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3 Gaussian Path Integrals and Determinants 



3.1 Preliminary Remarks 

The path integral expression for the propagator we have obtained is 

K(x f , Xi ; t f -t l )=N [ X(tf)=Xf D [x(t)]e (*/ft)S[x(t)] . (3.!) 

Jx(ti)=Xi 

We thus now need to make sense of, and develop rules for evaluating, such path integrals. 

For a general system described by an action S[x] an exact evaluation of the path integral 
is certainly too much to hope for. Indeed, even in the finite-dimensional case integrals 
of exponentials of elementary functions can typically be evaluated in closed form only in 
the purely quadratic (Gaussian, Fresnel) case, whereas more general integrals are then 
evaluated 'perturbatively' in terms of a generating function as in (A. 26). 

For path integrals, the situation is quite analogous. Typically, the path integrals that 
can be calculated in closed form are purely quadratic (Gaussian, Fresnel) integrals, i.e. 
actions of the general time-dependent harmonic oscillator type 

S[x] = ™£'(x(t) 2 -u,(t) 2 x(t) 2 ) , (3.2) 

perhaps with the addition of a 'source' term, 

S[x] ->• S[x] + J j{t)x(t) . (3.3) 

Then the strategy to deal with more general path integrals, corresponding e.g. to an 
action of the form 

S[x]= j\™x{tf-V{x{t)) , (3.4) 

is to reduce it to an expansion about a quadratic action. In practice this is achieved 
in one of two ways. Either the potential is of the harmonic oscillator form plus a 
perturbation, V(x) = Vq(x) + XW(x), and one defines the path integral via a power 
series expansion in A, a perturbative expansion. Or one defines the path integral by an 
expansion around a classical solution x c (t) of the equations of motion rax = —V'{x). 
To quadratic order in the 'quantum fluctuations' around the classical solution one then 
finds the action (3.2) with uj(t) 2 = ^V"(x c (t)). This turns out to lead to an expansion 
of the path integral in a power series in h, a semi- classical expansion. 

In either case, the Gaussian path integral can be evaluated in reasonably closed form 
and the complete path integral is then defined in terms of the generating functional 
associated with this quadratic action. In this section 3 we will deal exclusively with 
Gaussian integrals. The evaluation of more general integrals in terms of generating 
functionals will then be one of the subjects of section 4. 
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3.2 The Free Particle and the Normalisation Constant N 



We are now ready to tackle our first path integral. For obvious reasons we will consider 
the simplest dynamical system, namely the free particle, with Lagrangian 



L (x(t),±(t)) = -±(t) 2 . (3.5) 



We thus need to calculate 



K {x f , x t ; t f -U)=M f Xitf) _ Xf D[x(t)]e (WSo[x(t)] > (36) 

J x{ti) — Xi 

and we know that we should find the result (1.27), 



\2 



Since in this case we already know the result, we can use this calculation to determine 
the overall normalisation of the path integral from the continuum point of view (since 
M is universal: it depends only on m and (tf — U) and not on the potential V(x)). 

As described above, the general strategy is to expand the paths around a solution 
to the classical equations of motion. Here the starting action is already quadratic, 
but this expansion will have the added benefit of eliminating the boundary conditions 
x(tij) = Xij from the path integral. We thus split any path satisfying these boundary 
conditions into the sum of the classical path x c (t), 

x c (t) = x l + X ^^(t-t l ) , (3.8) 
and a quantum fluctuation y(t), 

x(t) = x c (t) + y(t) , (3.9) 
with y(t) satisfying the homogeneous boundary conditions 

y(t l )=y(t f )=0 . (3.10) 
Plugging this into the action, one finds 

S [x c + y} = S [x c ] + ^ /*' y(t) 2 (3.11) 
2 Jti 

where So[x c ] is the classical action, already given in (1.28), 

S [x c ] = m {x / ~ X 4 • (3-12) 
L J 2 (tf-U) V ' 
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There is no linear term in y(t) because we are expanding around a critical point of the 
action 5o[x(t)], and there are no higher than quadratic terms because the free particle 
action itself is quadratic. 

In the path integral, instead of integrating over all paths x(t) with the specified boundary 
conditions we now integrate over all paths y(t) with the boundary conditions (3.10). We 
thus find that the path integral expression for the kernel becomes 

K (x f ,x i ;t f -t i ) = e^ i /^ So ^M ( £>[y(t)]eHT Su ^ (3.13) 

Jy(ti)=y(t f )=0 

We thus automatically produce the classical action in the exponent, as in (3.7). 

To get a handle on the path integral over y(t), we integrate by parts in the action to 
obtain 



Jy(ti)=y{t f )=o Jy(ti)=y(tf)=o 

(3.14) 

Comparing with the fundamental Fresnel integral formula (A. 32), 



/ 

J — ' 



i d x JA a bX a x b /2 = ( i . _A_y 1/2 



d a xe ^ ah x x /z = (det^J , (3.15) 



we deduce that what this path integral formally calculates is the determinant of the 
differential operator (—df), 



[ D[y(t)]e^ £ y(t)(-9Mt) = ( Bet ™ [_^V 1/2 . (3.16) 

Jy(ti)=y(t f )=0 \ ZiriTi J 



I denote the (functional) determinant of an operator by Det to distinguish it from 
a standard finite-dimensional determinant. We will discuss this determinant in more 
detail momentarily. First of all, we note that comparison of (3.7) and (3.13) shows that 
the infinite normalisation constant N is related to this determinant by 

N = J J ^Det-^HflV 172 . (3.17) 

We also see that the normalised path integral measure D[y(t)] introduced in (2.28) is 
such that it normalises the free particle Gaussian fluctuation integral to unity, 

f D [y{t)] dfli , y(t)(-d?)y(t) = 1 (3 . 18) 

Jy(ti)=y(t f )=o 
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3.3 The Free Particle: Fluctuation Determinant and Mode Expansion 



Now let us return to the definition of the determinant of a differential operator. As in 
the finite-dimensional case, this determinant is formally defined as the product of the 
eigenvalues. That this is indeed what the path integral gives rise to can be seen more 
explicitly by expanding y(t) in normalised eigenmodes y n (t) of the operator (—df), 

y(t) = J2c n y n (t) , (3.19) 

n 

with 

(~ d t) y n (t) = x n y n (t) 

ftf 

/ yn(t)y m (t) — $m,n 
Jti 

Vn(U) = Vn(tf) = . (3.20) 

In terms of this decomposition, the action becomes 



m rt 



[ f y(t)(-d?)y(t) = ™Y, X nC 2 n ■ (3-21) 
Jti 2 



2 „, 

Thus the expansion of y(t) in terms of eigenmodes is tantamount to diagonalising the 
operator (this statement is obviously true more generally). Then the path integral over 
all y(t) becomes an integral over the c n , 

D[y(t)] = H( [ +0 °dc n ) , (3.22) 
iy(ti)=y(t f )=o n \J-oo J 

and thus the path integral reduces to an infinite product of finite-dimensional Gaussian 
integrals, with the result 

D[y(t)]eK2 J 

>y(ti)=y(t f )=o 



[ 

Jy( 



[ 

Jv( 



D[y (t)]^ ft! y{t){-d>)y{t) = U (f +0 ° dCn )MJ: n x n cl 

I „ \J—co / 



/ \ -1/2 

- (n^ A «) ■ < 3 - 23 > 

It is this infinite product which we have defined as the determinant. 

It will be useful for later to know explicitly the eigenvalues A n . The properly normalised 
eigenfunctions are 

I 2 t — U 

Vn(t) = 4 /- -sinnvr- l - . (3.24) 

y tf — ti tf — ti 

Since y- n (t) = —y n (t), the linearly independent solutions are y n (t) with n £ N and the 
corresponding eigenvalues are 

2 2 
TI TV 

K= w - ? . (3.25) 
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We thus have 



00 n 2 7r 2 



(3.26) 



This is clearly infinite, and thus Det _1 / 2 [— df] is zero, but this is compensated by the 
infinite normalisation constant in precisely such a way that one obtains the finite result 
(3-17). 

' " ' (3.27) 



V ~" ' 2irih l ~ l 'J V2mh(tf-U) 
For more comments on determinants and regularised determinants see sections 3.5 and 
3.6. 



3.4 The Harmonic Oscillator 



We have now accumulated all the techniques we need to tackle a more interesting ex- 
ample, namely the time-independent harmonic oscillator, with action 

S[x] = ^ f\x(t) 2 -^x(tf) . (3.28) 

Following the same strategy as for the free particle, we decompose the path into a 
classical path x c (t) and the fluctuation y(t), determine the classical action and the (still 
quadratic) action for y(t), and then perform the Gaussian integral over y(t). 

The classical path and action are (see Exercise 1) 



sin uj (tf -t) sin uj (t - U ) 



S[x c 



sin ujq (tf — ti 
m ujq 



sin uiq (tf — ti 



(xj + x 2 f) cosu;o(t/ — U) — 2xiXf 



2 sin wo(tf — t 
Expanding the action arund x c (t), one finds 

S[x c + y} = S[x c ] + ^ f f y(t)(-d 2 t - u 2 )y(t) 
2 Jti 

Thus the path integral we need to calculate is 

K{x f ,x i ;tf-t i ) = eW h W x - 



(3.29) 



(3.30) 



D[y(t)]enYlt f y(t)(-d 2 -^)y(t) 

y(ti)=y(tf)=o 

(3.31) 

This is once again a straightforward Gaussian (Fresnel) integral, and thus one finds, 
using the result (3.17), 



K(xf,xf,tf - U) 



I m Det 2%ft[-ffi Ji/h)S\Xr. 

2mh(tf-t l )\\Det^ E [-d 2 -u J 2 } 

Bet[-d 2 } (i/ti)S\x r ] 



m 



2mh(t f - U) V Bet[-df - uj 2 ] 



(3.32) 
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Since in writing the above we have taken into account the normalisation factor, the 
result should be well defined and finite. This is indeed the case. To calculate this ratio 
of determinants, we observe first of all that if the eigenvalues of the operator {—d 2 ) are 
X n , the eigenvalues fi n of the operator (— <9| — Uq) are \i n = \ n — Wq. Thus the ratio of 
determinants is 



Det[-a 



Det[-a t 2 -^ 2 ] 



v n A » 



1/2 



<4. 



-1/2 



(3.33) 



Even though this may not be obvious from this expression, the result is actually an 
elementary function. First of all, using the explicit expression for the X n , one has 



nc-?: 

n=l An 



n=l 



n 2 ir 2 



The function 



x 



n=l 



n 2 TT 2 ' 



(3.34) 



(3.35) 



is an even function of x with /(0) = 1 and simple zeros at x = ±nir. This shows (see 
also Appendix B) that this is an infinite product representation of 



/(*) 



sin a: 



Therefore, the final result for the ratio of determinants is 



Bet[-8 2 ] 



uo{tf- U) 



Y Det[— 8 2 — oj 2 ] y s'mujo(tf — U) 
and our final compact result for the propagator of the harmonic oscillator is 



(3.36) 



(3.37) 



K(xf,xf,tf - U) = 



m 



2irih(tf — ti) y sinwo(^/ — U 



Mtf-ti) (i/h)S[x r ] 



,(i/h)S[x c 



with 



S[x c 



2irih sin uq (tf — t 

(xf + x 2 f) cosoJo(tf — ti) — 2xiXf 



(3.38) 



(3.39) 



m coo 
2 sin wo (£/ — U) 

It is easy to see that this result reduces to that for the free particle in the limit uq — > 0. 
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Given the above result for the kernel, we also immediately obtain an expression for the 
partition function by setting X{ = Xf = x and integrating over x. This is a Fresnel 
integral, and the result is 

Z(t f - U) = ^— ^ . (3.40) 

2i sin 0( f 2 l > 

This can be expanded as 

oo 

Z ( tf _ t .j = J2 e -(i/h)(t f -ti)E n j (341) 

n=0 

where (cf. (2.33)) 

E n = hu Q {n+\) (3.42) 
are the harmonic oscillator energy levels. (Exercise) 



3.5 Gaussian Path Integrals and Determinants: the VVPM and GY For- 
mulae 

Frequently one encounters time-dependent harmonic oscillators with Hamiltonian 



m^l^^f** (3.43, 



and the quadratic action 

S[x} = ^j\x(t) 2 -u;(t) 2 x(t) 2 ) ■ (3-44) 

For instance, expanding the general action 

S[z] = % f f (x(t) 2 - *-V(x(t)) (3.45) 
^ Jti rn 

to second order around a classical solution x c {t) of the equations of motion mx = 
— V'(x), one finds the action (3.44) with 

co(t) 2 = ^V"(x&)) . (3.46) 

The resulting path integral is still Gaussian, and exactly the same strategy as above can 
be used to show that the path integral result for the kernel of the evolution operator is 
(cf. (3.32)) 



< x f \fe-W & dt 6 ^\x- >- I m J Det[ ~^ 2] e m)S[x c 

< fl 1 1 >_ V27Tih(t f -U)^Det[-d 2 -u(ty} 



(3.47) 
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Here x c {t) now denotes the classical harmonic oscillator solution with the given bound- 
ary condition, x c (ti) = x,- L and x c (tf) = Xf, S[x c ] is the classical action, and the fluctu- 
ation determinants are to be calculated for zero (Dirichlet) boundary conditions. 

In order to evaluate the result for the propagator (3.47), one needs to determine the 
classical action and the ratio of fluctuation determinants. The former is rather straight- 
forward provided that one can find the classical solution. An integration by parts shows 
that the classical action can be calculated in terms of the boundary values of x c (t) and 

S[x c ] = - jf (x c (t) 2 - u{tfx c {tf) = -[x f x c (t f ) - x % x c {ti)} . (3.48) 

The calculation of the ratio of determinants would be complicated if one tried to cal- 
culate these determinants directly, as we did in the time-independent case. Fortunately 
there are two elegant shortcuts to calculating this ratio of determinants which are finite- 
dimensional in nature and do not require the calculation of a fuctional determinant. I 
will briefly describe these below. 

One can for instance use the useful and remarkable result that the ratio of determinants 
can be calculated from the classical action via the Van Vleck - Pauli - Morette (VVPM) 
formula 

Det[-a t 2 ] _ i / d*s[x c ] (3 _ 4g) 



m 



V 2tu/% - U) V Det[-d t 2 - oj(t) 2 ] y^ih V d Xi dx f 
More generally, for a d-dimensional quantum system, the 2nd derivative of the classical 
action would be replaced by the d-dimensional VVPM determinant 

'd 2 s[x c 



^ ^det 



dxfdx^ 



dxidxf 

This is a non-trivial but standard and well-known result. Notice that, to evaluate the 
ratio of quantum fluctuations in this manner, one only needs to know the classical 
action. 

In the case of the harmonic oscillator with constant frequency, agreement between the 
VVPM formula and the result we obtained in (3.38) can be immediately verified from 
the classical action (3.39) which gives 

d 2 S[x c ] mujQ 



dxidxf sin (t/ — ti) 



(3.50) 



Alternatively, instead of the VVPM result one can use the equally remarkable, but 
apparently much less well known, Gelfand-Yaglom (GY) formula that states that 

Deth^-c^) 2 ] _ FUtf) 

Det[-df] " tf-U 1 ' 
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where F^it) is the solution of the classical harmonic oscillator equation 

(d*+u(t) 2 )F u (t) = (3.52) 

with the initial conditions 

F u (U) = F w (ti) = l . (3.53) 
Thus one has the simple result 



/ m / Det[— df] _ I m 

2mU{t f - U) V Det[-<9 t 2 - uj{t) 2 ] ~ V 2mhF u (t f ) 



(3.54) 



It is quite remarkable that the ratio of fluctuation determinants, which involves the 
product over all eigenvalues of the operator — {d 2 + ui(t) 2 ) (with zero boundary condi- 
tions), can be expressed in terms of the zero mode (solution with zero eigenvalue) of 
the same operator with the GY boundary conditions (3.53). 

Once again, as an example we consider the constant frequency harmonic oscillator. The 
solution of the classical equations of motion satisfying the GY boundary conditions is 
evidently 

F^t) = — sin w (t - U) ■ (3-55) 

Thus 

F Wo (tf) = — smoj (tf - U) (3.56) 



and the GY formula predicts 

Det[— d\ — uj 2 ] _ sin (t/ — U 



(3.57) 



Det[-9 t 2 ] u (t f -ti) ' 

in perfect agreement with the result (3.37). 

Some comments: 

1. If f(t) is a solution of the harmonic oscillator equation with f(U) = 0, the GY 
solution is evidently F w (t) = f(t)/f(U) with 

F ^f) = 777T • ( 3 - 58 ) 

If, on the other hand, f(t) is any solution of the harmonic oscillator equation with 
f(U)^0, the GY solution can be constructed from it as 

F ul (t) = f(U)f(t)£dt'j^ (3.59) 

In particular, one has 

F u (tf) = f(U)f(tf) £ dt . (3.60) 
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2. The denominator tf — U in the GY formula (3.51) can be interpreted as Fo(i/), 
where 

F (t) = t-ti (3.61) 

solves the free-particle (lo = 0) equations of motion with the GY boundary condi- 
tions, so that one can write (3.51) in the more suggestive form 

Det[-d t 2 - cj(t) 2 } _ F u (t f ) 



Bet[-d 2 } F (tf) ■ 

3. Comparison with the VVPM formula gives the relation 

d 2 S[x c ] _ m 



(3.62) 



(3.63) 



dxidxf F u (tf) 

With a bit of effort this purely classical (in the sense of classical mechanics) formula 
can be proved directly via Hamilton- Jacobi theory (see Appendix C), but this by 
itself provides no insight into the reason for the validity of either the VVPM or 
the GY formula. 

A slick proof of the GY formula, in the form (3.62), has been given by S. Coleman in 
his Erice lectures on The uses of instantons (Appendix A), reprinted in 

S. Coleman, Aspects of Symmetry, Selected Erice Lectures, Cambridge University Press 
(1985). 

This proof is also reproduced in 

L.S. Schulman, Techniques and Applications of Path Integration, Dover (2005). 

It works roughly as follows (to make this argument more precise one should insert words 
like "Fredholm operators" etc. in appropriate places): 

Let Fu t \(t) be the solution of the equation 

(-d 2 - u(t) 2 )F^ x (t) = XF u , x (t) (3.64) 
with the GY boundary conditions 

F U , X (U) = , F UtX (U) = 1 • (3-65) 

Thus F w fi(t) is what we called F u (t) above, and F u \(t) has the property that F u \(tf) = 
iff A is an eigenvalue of the operator (—d 2 — oj{t) 2 ) with zero Dirichlet boundary 
conditions (because then F Ut \(i) is the corresponding eigenfunction) . 

The claim is now that 

Det[-ff-q,(t) 2 -A] _ F u>x (t f ) 

Bet[-d 2 - A] F Mt f ) ^ 
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for all A € C. Here, Det is again defined to be the product of all eigenvalues. In 
particular, this implies the GY result (3.62) for A = 0. 

To establish this claim, one considers the left and right hand sides as functions of the 
complex variable A. The left hand side is a meromorphic function of A with a simple 
zero at each eigenvalue X n of (— d 2 — u:(t) 2 ) (an eigenvalue \ n of (— d 2 — w(i) 2 ) is a zero 
eigenvalue of (—d 2 — u>(t) 2 ) — A n )), and a simple pole at each eigenvalue Ao, n of (—8%) 
(for the same reason). By the remark above (3.66), exactly the same is true of the right 
hand side. In particular, the ratio of the left and the right hand sides has no poles and 
is therefore an analytic function of A. 

Moreover, provided that oj(t) is a bounded function oft, for A sufficiently large, |A| — > oo, 
one can ignore oj(t), and hence both the left and the right hand side go to 1 in that limit 
(everywhere except on the real positive line where one can find large real eigenvalues). 

Putting these two observations together, one concludes that the ratio of the two sides 
is an analytic function of A that goes to 1 in any direction except perhaps along the 
positive real axis, and this implies that the ratio is equal to 1 identically, which concludes 
the proof of the identity (3.66). 

Another elegant continuum (i.e. non-discretised) proof of the GY formula (3.51) can be 
found in 

H. Kleinert, A. Chervyakov, Simple Explicit Formulas for Gaussian Path Integrals with 
Time- Dependent Frequencies, Phys. Lett. A245 (1998) 345-357; quant -ph/9803016. 

Yet another proof can be assembled from sections 3.3, 3.5 and 34.2 of 

J. Zinn-Justin, Quantum Field Theory and Critical Phenomena (Oxford Science Publi- 
cations, 1989). 

A proof along these lines is also outlined in Appendix A of 

R. Rajaraman, Solitons and Instantons (North Holland 1982, 1987) 

3.6 Some Comments on the Regularisation of Determinants 

In the above we have repeatedly encountered determinants of infinite-dimensional oper- 
ators, and we have treated them in quite a formal and cavalier way. In the following, I 
want to (very) briefly indicate some ways to define these objects in a more satisfactory 
manner. 

There are certain infinite-dimensional operators for which the definition of a determinant 
poses no real problem. For example, for most intents and purposes, trace class operators 
K (i.e. operators for which the trace exists), and operators of the form I + K, where I is 
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the identity operator and K a trace class operator, behave like finite-dimensional linear 
operators (matrices). For invertible (n x n)-matrices M one can write the determinant 
(with K = I - M) as 

det M = det(I + K) = 1 + tr K + tr( A 2 K) + ... + tv(A n K) . (3.67) 

Here A P K (the p'th anti-symmetric power of K) denotes the operator implementing 
the induced action of K on anti-symmetric p-vectors. Analogously, one can define the 
(Fredholm) determinant of an invertible operator I + K with K trace class by 

Det F {I + K) = 1 + TrK + Tt(A 2 K) + ... , (3.68) 

since this series is absolutely convergent. However, most operators appearing in physics 
are not of this form, and hence one needs to be more creative. 

In section 3.4 we had seen that, even though Det(— d 2 — coq) (the product of the eigen- 
values) diverges, the ratio of this determinant and the (equally divergent) free particle 
determinant gave us a finite result. More generally, the ratio of determinants that ap- 
pears in (3.47) can be interpreted as defining a regularised functional determinant of 
the operator (— d 2 — u(t) 2 ), in the sense of 

^ r o2 / n2i Det[-a 2 - Uj(t) 2 } 

Bet reg [-d 2 - u>{t) 2 ] := Det[ _ffi (3-69) 

It is this regularised determinant that the normalised path integral with measure D[x] 
computes (cf. (2.28,3.18)), 

^(t)]e^A f ym-^-^)y(t) = ( DeW _ 9 2 _ u{t?] y^ 2 . 

y(ti)=y(t f )=0 V ! 

(3.70) 

However, usually in the physics literature one adopts a slighly different attitude. Instead 
of regularising explicitly by means of the free particle determinant (which is neverthe- 
less natural from the path integral point of view) one attempts to define meaningful 
individual (instead of ratios of) regularised functional determinants in other ways, e.g. 
via the so-called £-function or heat kernel regularisation. 

For a sufficiently reasonable (elliptic, self-adjoint, . . . ) operator A with eigenvalues A n , 
define the spectral (-function of A to be the function 

U(s) = J2K S ■ (3-71) 

n 

This generalises the definition of the ordinary Riemann ("-function ((s), 

oo 

C(s) = £ n- . (3.72) 

n=l 
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To see the relation between the spectral function and the determinant, one differenti- 
ates once, 

d(*) = -E n ~ sl °g A « > ( 3 - 73 ) 

n 

to conclude that 

Ca(0) = - E lQ g X - = ~ log LI X n ■ (3-74) 

n n 

Formally, therefore, one has 

nA„ = e-Cl(0) (3J5) 

n 

but, at this point, the left hand side is as ill-defined as the right hand side because 
Ca(s), as it stands, will be convergent only for Re(s) sufficiently large and positive (and 
not at s = 0). 

Likewise, the ordinary Riemann (^-function, as it stands, converges only for Re(s) > 1. 
However, in that case it is well known that ((s) can be analytically continued to a 
meromorphic function of s in the entire complex s-plane, with a pole only at s = 1. In 
particular, one then has cute results like 

C(0) — | (£~ii) 

C(-l) = -£ ( En=in ) (3.76) 
C("2) = (En=in 2 ) 

(for a poor man's (handwaving) proof of these identities, see Appendix D) as well as 
(we will need this below) 

C / (0) = -ilog2vr . (3.77) 

Note that, with 

oo 

C'(0) = -X>g" , (3-78) 

71=1 

this result can be written as the equally cute and astounding identity 

oo 

Y[ n 2 = 2vr . (3.79) 

n=l 

Analogously, under favourable circumstances the spectral £- function (a(s) can be ex- 
tended to a meromorphic function of s which is holomorphic at s = 0, and then (3.75) 
can be used to define the C-function regularised determinant of A via 

Det c A = e~&(°) . (3.80) 

To illustrate this method, let us go back to the calculation of the free particle determi- 
nant in section 3.3. There we had seen that the eigenvalues of the operator A = —df 
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(with Dirichlet boundary conditions) are 

K = ^ r n=l,2,... (3.81) 
where T = tf — U. To define the determinant, we construct the spectral (- function 

C*M = t (f)'" = © 2 "« 2 »> = ' 2slOS{T/ ' ) C(^ (3.82) 

and calculate 

C' A (0) = 2C(0) log(T/vr) + 2C'(0) = - log(T/vr) - log(27r) = - log(2T) . (3.83) 
Thus we conclude that, with ("-function regularisation, 

Det c [-df] = 2T . (3.84) 
Likewise, for the harmonic oscillator determinant of section 3.4, with 

2 2 

A„ = ^--w 2 n = l,2,... (3.85) 



one finds 



In particular, we have 



Be k [-df-^]=2T S ^^ . (3.86) 



Det^[— dl — ujq] sin Tw 



Det c [-d t 2 ] Tw 
in perfect agreement with the previously obtained result (3.37), 



Det[-d 2 ] Mtf-U) 



Det[— df — ujq] y sin wo (*/ — *i) 



(3.87) 



(3.88) 



While this (-function regularised definition of the determinant captures most of the 
essential properties of the standard determinant, some care is required in manipulating 
these objects. For example, one important property of the standard determinant in 
the finite-dimensional case is its multiplicativity det(MiV) = det(M) det(iV). To what 
extent an analogous identity 

Det c (AB) = Detf A Det f B (3.89) 
holds in the infinite-dimensional case is discussed e.g. in elemntary terms in 

hep-th/9804118; E. Elizalde, On the concept of determinant for the differential opera- 
tors of Quantum Physics, hep-th/9906229, 

and, from a more mathematical point of view, in (warning: not for the faint of heart) 

M. Kontsevich, S. Vishik, Determinants of elliptic pseudo- differential operators (155 p., 
hep-th/9404046). 
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3.7 Exercises 



1. Verify (3.29). Instead of calculating the classical action by integration, try to 
determine it in an intelligent way by proving first that only boundary terms con- 
tribute to the classical action, 



2. Verify that the harmonic oscillator kernel (3.38) 

• satisfies the harmonic oscillator Schrodinger equation (1.15) 

• and the initial condition (1.14), 

• and reduces to the free particle propagator for ujq — > 0. 

3. Verify the result (3.40) for the partition function of the harmonic oscillator, and 
the expansion (3.41). 

4. Fill in the missing steps in the proof of the identity (3.63) given in Appendix C. 
In particular, verify (C.10). 

5. Find an elementary proof of the VVPM formula 



S [ x c] = -j[XfX c (tf) ~ XiX c (U)] . 



(3.90) 



Det[-d 2 ] 




(3.91) 



Dct[-a t 2 -cv(t) 2 ] 



or the GY formula 



Det[-9 t 2 ] 



Fo(tf) 



(3.92) 



Det[-9 t 2 -uj{t) 2 } 
and publish it and/or tell me about it. 
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4 Generating Functionals and Perturbative Expansions 



In this section we will study some other important properties of the path integral, in 
particular the perturbative and semi-classical expansion of non-Gaussian path integrals. 
The treatment in this section is somewhat more cursory than in the previous sections - 
the main intention is to give a flavour of things to come (in a course on quantum field 
theory, say). 



4.1 The Generating Functional Z\J\ 

The main objects of interest in quantum field theory are vacuum expectation values 
of time-ordered products of field operators. These matrix elements can be obtained 
from a generating functional which, in turn, can be expressed as a path integral. This 
motivates the following discussion of these concepts in the quantum mechanical context. 

Before turning to the path integral, we introduce the generating functionals for the 
correlation functions (n-point functions) 

Gfi(x(ti), . . . ,x(t n )) = <Xf,tf\T(x(ti)...x(t n ))\xi,ti> 
G(x(ti), . . . ,x(t n )) = <0\T(x(t 1 )...x(t n ))\0> . (4.1) 

They are defined as 

(i/h) n r*f 



(in) rf 

z fi[j] = E^tW dt 1 ...dt n j(t 1 )...j(t n )G fi {x(t 1 ),...,x(t n )) 

n=0 n - Jt ' 

= <x f ,t f \Te^ n ^U dtj(t)x(t) lx ^ ti> (42) 
Z \j\ = E^lM dt 1 ...dt n j(t 1 )...j(t n )G(x(t 1 ),...,x(t n )) 

n=0 U - J ~°° 

= <0\re^/ h ^-^ dtj ^ £ ^\0> . (4.3) 



From these generating functionals, the individual correlation functions can evidently be 
reconstructed by differentiation, 

« - (D' J'Su ^ 

G Wtl) ,.. .,*(*„)> = g)' w g [g bH • («> 

For Zfi[j] we can easily deduce a path intgral representation. Using (2.60) and the 
definition (4.2), one finds 

Z fi [j] = N f X{tf)=Xf D[x(t)\ e (i/K)S[x(t);j(t); t f ,U] ) (45) 
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where t 

S[x{ty,j(t);t f ,t i ] = S[x(t);t f ,t i }+ I '* dt j(t)x{t) (4.6) 

Ju 

is the action with a source term (or the action with a coupling of x(t) to the current 
j(t)). Thus the generating functional Zfi[j] is the path integral for the action with a 
source term. 

Our aim is now to find a similar path integral representation for Z[j]. For that we need 
to project from the states \x,t > to the ground state |0 >. To that end we first expand 
the state \x, t > in a basis of eigenstates \n > of the Hamiltonian, 

|x,t>=e(W^|x>=^e(W^|n><n|x>=^e(W^<(x)|n> . (4.7) 

n n 

Thus for the correlation functions (2.60) of time-ordered products of operators one finds 
Gfi(x(ti), . . . ,x(t n )) = ^e-W^/^ + ^^x 

m,n 

^n{x f W m {xi) <n\T{x{ti)...x{t n ))\m> . (4.8) 

To accomplish the projection onto the vacuum expectation value, we now take the limit 
tf^ — > ±oo. This can be understood in a number of related ways. They all amount 
to the statement that, in the sense of distributions, exp(-itE) — > for t — > oo, the 
dominant contribution in that limit coming from the smallest possible value of E, i.e. 
from the ground state. Explicitly, one can for instance replace E n — > (1 — ie)E n for 
a small positive e and then take the limit e — > at the end. Alternatively, one can 
"analytically continue" to imaginary time, take the limit there, and then continue back 
to real time. In whichever way one proceeds, one can conclude that 

G{x(h), . . . , x(t n )) = < 0\T(x(h) . . . x(t n ))\0 > (4.9) 

(i/h)(tf - U)E r x(t f )=x f /-/%\a\ u\ + +1 

= I™ , , \— : / D[x(t)]x(t 1 )...x(t n )e( t /^ S i x (^ t f^ . 

As the left hand side is independent of the boundary conditions x imposed at t — > ±oo, 
so is the right hand side. Passing now to the generating functional Z[j] (4.3), we can 
once again rewrite the infinite sum as an exponential in the path integral to deduce 
(suppressing the dependence on the boundary conditions) that 

Z[j] ~ J D[x(t)} e {i/K)S[x(t);j(t); t u = ±oo] j (4 10) 

where 

/+oo 
dt[L(x(t),x(t))+j(t)x(t)\ . (4.11) 
-oo 
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The proportionality factor (normalisation constant) is fixed by 

Z[j = 0] =< 0|0 >= 1 . (4.12) 

We therefore obtain the final result 

/D[x(t)]e(^)^);^);^ = ±oo] 
m ~ lD[x{t)]e m)S[x(t);t u = ±oo] ■ { ^ l6) 

The right hand side is independent of the boundary conditions imposed at t = ±oo 
provided that one chooses the same boundary conditions in the numerator and the 
denominator. 



4.2 Green's Functions and the Generating Functional for Quadratic 
Theories 



As we will reduce the calculation of a general path integral and its generating functional 
to that of a quadratic theory, in this section we will determine explicitly the generating 
functional for the latter. 

For finite-dimensional Fresnel integrals one has (see e.g. Exercise 1.5.3 and equation 
(A.13)) 

JA ab x a x b /2 + ij a x a = ( det A\ V2 fi -iG ab j a j b /2 



J d d x e' 1 



2m J 



(4.14) 



where G ab is the inverse matrix ("Green's function") to A ab , G ab A bc = 5 a c . Thus the 
"generating function" is 



J d^X Q^^abX a X b /2 + ij a x a 



= e 

It follows from this that the "2-point function" 



J d^X Q^-^abX a X /2 

iG ab j a j b /2 



< x c x d >: 



/ d d x x c x d e iA a bX a x /2 

j d^X Q^^a,b xCl X b 1 2 



IS 



< x c x d > 



z o[j\ 



1 d 1 d 

i dj c i dj d 

iG cd . 



3=0 



3=0 



(4.15) 



(4.16) 



(4.17) 
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Thus the "2-point function" is the "Green's function" of the Fresnel integral. 

Higher moments (n-point functions) can be calculated in a similar way. For n odd they 
are manifestly zero. For the 4-point function one finds (Exercise) 

< X a XbX c Xd > = < X a Xb X x c x d > + < x a x c X x b x d > + < X a Xd X x b x c > (4.18) 

etc. The general result, expressing the 2n-point functions as a sum over all possible 
pairings P(x 1 , . . .,x 2n ), 

<X 1 ...X 2n >= < x h x 'i2 > ■ ■ ■ < x i 2n -i x 'i2n > ( 4 - 19 ) 

P(xi,...,X 2n ) 

(there are (2n — 1)!! terms) is also easily deduced from the generating function zo[j]. 
In the quantum field theory context, this result is known as Wick's Theorem and, even 
though a simple result, is of enormous practical significance in perturbative calculations. 

We now consider the analagous question for harmonic oscillator path integrals. In this 
case, one finds, in precise analogy with the finite-dimensional case, 

where G(t,t') is a Green's function (the inverse) of the operator —d$ — w(i) 2 , 

(-df-co(tf)G (t,t') = S(t-t') . (4.21) 

For the finite-time path integral, the Green's function that appears here would have 
been determined by the Dirichlet (zero) boundary conditions at tij to be the Green's 
function with 

G o (t f ,t') = Go(t,ti) = . (4.22) 

In the present case (infinite time interval), the relevant Green's function is implicitly 
determined by an ie prescription. In particular, for the time-independent harmonic 
oscillator with constant frequency ujq one finds 

Go(M') = -Le-^ol* - *'! . (4.23) 

Once we know the generating functional, we can use it to calculate n-point functions. 
In particular, for the two-point function one has 

'h 5 h 5 

7iJJ{t)~iWt' 
ih 

= — ( 

m 

That the time-ordered product < 0\T(x(t)x(t'))\0 > is a Green's function can of course 
also be verified directly in the standard operator formulation of quantum mechanics 
(Exercise). Likewise, higher n-point functions can be expressed in terms of products of 
2-point functions (Wick's theorem again). 



< 0|T(x(*)z(O)|0 > 



J i=o 

-Go(M') • (4.24) 
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4.3 Perturbative Expansion and Generating Functionals 

As we have seen, the generating functional Z[j] encodes all the information about the 
n-point functions G(xi, . . . ,x n ). However, this is only useful if supplemented by a 
prescription for how to calculate Z[j]. Since in practice the only path integrals that 
we can do explicitly are Gaussian path integrals and their close relatives, the question 
arises how to reduce the evaluation of Z[j] to an evaluation of Gaussian integrals. 
This is achieved via a perturbative expansion of the path integral around a quadratic 
(Gaussian) action. The (assumed) small perturbation expansion parameter can be a 
coupling constant A, as in 



with Hq a harmonic oscillator Hamiltonian. In this case one is studying a path integral 
counterpart of standard quantum mechanical perturbation theory. 

Alternatively, the small parameter could be Planck's constant Ti itself. In this case one 
is interested in evaluating the path integral 



for Ti — > as a power series in % by expanding the action around a classical solution x c (t). 
This is a semi-classical expansion of the path integral, a counterpart of the standard 
WKB approximation of quantum mechanics. 

Both physically and technically (in one case one has a small parameter in front of a part 
of the action, the perturbation, in the other a large parameter l/fi in front of the entire 
action) these two expansions appear to be quite distinct. Calculationally, however, 
they are rather similar, since in both cases the path integral can be reduced to a series 
expansion in derivatives of the generating functional of the quadratic theory. This is 
immediate for the perturbative A-expansion (which we will consider in this section) 
but requires a minor bit of trickery for the ^-expansion (hence the excursion into the 
stationary phase approximation for finite-dimensional integrals in section 4.5). 

For definiteness, we will assume that the perturbation W arises from a velocity-independent 
perturbation of the potential 



with Vq{x) a harmonic oscillator potential. For more complicated perturbations one 
would have to go back to the phase space path integral, introduce sources for both x(t) 
and p(t), etc. Then the action takes the form 



H = H + \W , 



(4.25) 




(4.26) 



V{x) = V (x) + XW(x) , 



(4.27) 




(4.28) 
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where Sq[x] is the free action, and Si[x] the perturbation or interaction term. 

To calculate the path integral, one introduces a source-term for the free action Sq[x], 



So[x,j] = S [x] + j dt j(t)x(t); , (4.29) 



and determines ^/i,o[j] or Zo\j]. Focussing on the latter, the vacuum generating func- 
tional Z[j] for the perturbed action can then be written in terms of that for the free 
action as 

Z[j]=Me {iX/h)Sl[ ^ ] Z [j] . (4.30) 
The normalisation constant, determined by the condition Z[j = 0] = 1, is 



AT 1 



Zo[j] 



(4.31) 



3=0 



This result for Z[j] is manifestly a power series expansion in A - and the fact that this 
perturbative expansion is so straightforward to obtain in the path integral formalism is 
one of the reasons that makes the path integral approach to quantum field theory so 
powerful. 

Moreover, using the explicit expression (4.20) for the generating functional Zq[j\ in terms 
of the Green's functions of the free theory, one sees that the generating functional Z[j], 
and thus all the n-point functions of the perturbed theory, are expressed as a series 
expansion in terms of the Green's functions of the unperturbed theory. A graphical 
representation of this expansion leads to the Feynman diagram expansion of quantum 
mechanics (and quantum field theory). 



4.4 The Stationary Phase Approximation for Oscillatory Integrals 

The integrals of interest in this section are oscillatory integrals of the kind 

f + °° dxeWWW . (4.32) 

The basic tenet of the stationary phase approximation of such integrals is that for small 
h, h — >■ 0, the integrand oscillates so rapidly that the integral over any small x-interval 
will give zero unless one is close to a critical point x c of f(x), f'(x c ) = 0, for which to 
first order around x c there are no oscillations. This suggests that for h — > the integral 
is dominated by the contribution from the neighbourhood of some critical point (s) x c 
of f(x), and that therefore in this limit the dominant contribution to the integral can 
be obtained by a Taylor expansion of f(x) around x = x c . 
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To set the stage for this discussion, we will first reconsider the Gaussian and Fresnel 
integrals of Appendix A from this point of view. There we had obtained the formula 
(A.13), 



/ 



°° dx e^/2 + ijx = Pjl e ~if/2a (4 . 3 3) 



by analytic continuation from the corresponding Gaussian integral and/or completing 
the square. A (for the following) more instructive way of obtaining this result is to 
consider the integral 

'' dx eW h M x ) (4.34) 



— oo 

2 , 



for some quadratic function of x (q(x) = ax 2 /2 + jx + c, say). Such a function has a 
unique critical point x c (x c = —j/a), and since q(x) is quadratic and q'(x c ) = one can 
write q(x) as 

q{x) = q(x c ) + \{x-x c ) 2 q"{x c ) . (4.35) 

Thus the integral is 



+ " m)q(x) = Ji/h)q(x c ) [ + °° Hllp m)q»(x c )yy2 



I ( i,. v \'i i>><i\,n ... l A'/")<n-'-ri j 

J —oo 



C 



(i/h)q(x c ) / M_ (436) 



This reproduces and generalises (4.33). 

It will be more convenient to move the factor V^iriTi to the left, so that the result is 

_^L_ f + °° drr p (i/h)q(x) = 1 p (t/n)g(a c ) _ (437) 
V^TrihJ-oo yq"{ x c) 

We see that a Fresnel integral is determined exactly by the contribution from the critical 
point and the quadratic fluctuations around it. Another way of saying this is that, as 
we will see, for a Gaussian/Fresnel integral the stationary phase approximation is exact. 

Now let us consider the oscillatory integral 

/[/] = -=L= / + °° dx e (*/&)/(*) (4.38) 

for a general function f(x). Let x c denote the critical point of f(x) with the smallest 
absolute value, and assume that x c is isolated. Then we can expand f(x) around x c as 

f{x) = f{x c ) + \{x-x c ) 2 f'\x c )+R(x-x c ) , (4.39) 

where the rest R(x — x c ) is at least cubic in (x — x c ). Thus the integral is 

/[/]= 1 e (i/h)f(x c ) f + °° dxe (i/h)[f"(x c )(x-x c f/2 + R(x-x c )} (44Q) 
y/2mh J-oo 
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The stationary phase approximation (also known as the saddle point approximation) to 
the integral amounts to ignoring the higher than quadratic terms encoded in R(x — x 
and leads to the approximate result 



f 00 dx e mm*) „ . 1 e mm*) . {4A1) 

J-00 ^fWixT) V ' 



_ ; +00 dxe(^)/(-)._i 

To justify this approximation, one needs to show that the contributions due to the 
remainder R{x — x c ) are indeed subleading in % as K — > 0. We will establish this below. 

Another cause for concern may be that, while we have (hand-wavingly) argued that the 
dominant contribution to the oscillatory integral should arise from a small neighbour- 
hood of the critical point(s), in order to arrive at (4.41) we have taken the integral not 
over a small neighbourhood of x c but, quite on the contrary, over (—00, +00). 

To justify this, consider the contribution to the integral from an interval [a, b] without 
critical points. In that interval, one can change the integration variable from x to f{x) 
(this would not be allowed if f(x) had a critical point in the interval). Then it is easy 
to see (e.g. by an integration by parts) that the integral 

r b 



I 

J a 



dxe (i/h)f(x) = ( h j _ ( 442 ) 



Thus regions without critical points contribute 0(K) terms to the integral. On the 
other hand, (4.41) shows that, regardless of what the neglected terms do, there are 
some contributions to the total integral which are of order 0(Ti 1 / 2 ). Thus these must 
be due to the contributions from (integrals over arbitrarily small neighbourhoods of) 
the critical points. As these are dominant relative to the 0(K) contributions as % — > 0, 
the difference betwen integrating over such a neighbourhood of the critical point and 
integrating over all x is negligible in this limit. 

To analyse the contributions due to R(x — x c ) and to make the dependence of the various 
terms on Ti more transparent, it is convenient to define the fluctuation variable y not as 
(x — x c ), as was implicitly done in (4.36), but via 

x = x c + Vhy . (4.43) 

This has the effect of making the Gaussian part of the integral independent of h, 

(x-x c ) 2 /2fi = y 2 /2 . (4.44) 

Moreover, since R(x — x c ) is at least cubic, (\/h)R{yfr,y) is now a power series in strictly 
positive powers of Vh, 



r(y) = (l/h)R(Vhy) = h 1 ' 2 (x c )y 3 /3! + hf^(x c )y 4 /4\ + ... (4.45) 
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even (odd) powers of y appearing with integral (half-integral) powers of H. All in all, 
we thus have 

1 f °° dx e WWW = -J_ e (i/ft)/(*c) [ + °° dy jW'iPcW ! /2 + r(y)] _ 

(4.46) 



To obtain the higher order corrections to the stationary phase approximation, one can 
expand expir(y). Remembering that only even powers of y in that expansion give a 
non-zero contribution to the y-integral, one concludes that this expresses the integral 
as a power series in (integral) powers of h, 

-L= [ + °° dx e = - % e WWfrc) (l + . .) + h\. . .) + . . .) (4.47) 

yl-KlkJ-oo Vf"( x c) 

Alternatively, this series can be expressed in terms of generating functions as 

1 + . .) + h 2 (. ..) + ... = e ir (-^/^)e-^ 2 /2/"(^)| j=0 . (4.48) 



The stationary phase approximation can be used to calculate the integral with reason- 
able accuracy for very small h, U — > 0. However, the series expansion should not be 
expected to converge in general, and the series is only an asymptotic series. If contribu- 
tions from all critical points are included, it is under certain conditions possible to obtain 
error estimates, but in practice applications of the stationary phase approximation are 
usually restricted to (4.41). 



4.5 The Semi-Classical Approximation 



It is now straightforward to formally apply this stationary phase approximation to 
path integrals. The crucial point is that in the path integral formulation of quantum 
mechanics the classical limit and the emergence of classical mechanics from quantum 
mechanics are extremely transparent: in the limit Ti — > 0, the path integral is dominated 
by paths that are critical points of the action, i.e. the classical solutions to the equations 
of motion. This should be contrasted with the much more cumbersome eikonal or WKB 
semi-classical approximation to the Schrodinger equation. 

Essentially all the work has already been done in sections 3.5 and 4.6 and we can be 
brief about this here. Instead of the function f(x) we have an action S[x], and the 
semi-clasical approximation amounts to expanding S[x] around a critical point, i.e. a 
classical solution x c (t) of the corresponding Euler-Lagrange equations. This is precisely 
the procedure we had already advocated for how to deal with general non-Gaussian 
path integrals. We now see that this will lead to an expansion of the path integral in 
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powers of h, the stationary phase approximation to the path integral agreeing with the 
harmonic oscillator path integral of section 3.5. 

First of all, with 

x (t) = Xc (t) + Sx(t) , (4.49) 

we expand the action as 

r*/ , 5S 



rtf fig 

S[x] = S[x c ] + ^ dt ^—.\ x{t)=Xc{t) Sx(t) 

1 /■*/ , /■*/ 5 2 S 
+- / dt dt' - , xr . 

2 7 ti 4 Sx(t)Sx(t 



= S[x c ] + -JJ dt JJ dt' g^jfa^jTj \x(t)=x c (t)^ x (t)^x(t') + ... (4.50) 



For an action of the standard form 



S[x} = ^ [ tf (x(tf-±V(x(t)) (4.51) 



the quadratic term is 



2 * ^ l*w-c(t)feWfe(0 = f / V(*) 2 - ^V"(x c (t))fe(t) 2 ) 

(4-52) 

As already noted in section 3.5, this is the action of a harmonic oscillator with time- 
dependent frequency 

co(t) 2 = ^V"{x c {t)) . (4.53) 

Thus the stationary phase or semi-classical approximation to the path integral is (cf. 
(3.47)) 



r (tf) ~ Xf Dixit)] ei^i^f^ « / m Det[-ff] {i /K)S[x c ] 

J x (U)= Xi pmHtf-t^Detl-dt-uitrf 

(4.54) 

This can be evaluated using either the VVPM or the GY method. 

The difference between (3.47)) and (4.54) is that in the former case one was dealing 
with a quadratic action and the result was exact (i.e. the semi-classical approximation 
is exact for the harmonic oscillator), whereas here this is really just the semi-classical 
approximation. Note also that x c in (4.54) refers to a classical solution of the full 
equations of motion mx = —V'(x) whereas x c in (3.47) is of course a solution of the 
harmonic oscillator equation. 

If desired, higher order corrections in U to the semi-classical result can be calculated, as 
in (4.48), by using the generating functional of the quadratic theory. 
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4.6 Scattering Theory and the Path Integral 

Formal scattering theory can be succinctly described by the evolution operator Uj(tf,ti) 
in the interaction representation. Thus let the Hamiltonian be (in the simplest case of 
two-body potential scattering) 

H = H + V (4.55) 
where Hq is the free-particle Hamiltonian Ho = p 2 /2m. Then Ui(tf,ti) is given by 

Ui(t f ,U) =e(W/#o e -(i/h)(t f -U)H e -{i/h)tiH (456) 

The M0ller operators 

U± = lim Z7/(0, t) (4.57) 

t->-=Foo 

act on free particle states \k a > (plane wave eigenstates of the free Hamiltonian Ho) 

<x\k a >=j^ W ~ 2 e i ^ , (4.58) 

p a = hk a , Ho\k a >= E a \k a >, to produce the stationary scattering states Y>J(x); 

n±\k a >=\^> . (4.59) 

The S-matrix elements S ab are the transition amplitudes among the asymptotic in and 
out scattering states >, 

S ab =<^t > , (4.60) 
and therefore they are the matrix elements 

S ab =< ka\S\h > (4.61) 

of the scattering operator 

S = riLn + = lim C/>(*/,ti) (4.62) 

t f^—t^OG 

between free particle states. 

Such matrix elements can readily be expressed in terms of the path integral. First of 
all, we have 

S ab = lim ^{i/^tfHo e -{i/n){t f -U)H e -{i/K)UHo^ b> 

tf,i— >±oo 

= lim e (i/h)(t f E a -UE b ) < ^,-(i/n)(t f -t i )H^ b> (463) 

t/,i-*-±oo 
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To pass from the momentum space matrix elements of the evolution operator for H to 
the kernel (the position space matrix elements), we perform a double Fourier transform 
(this is a special case of (2.61)), 

< k a \e-( i /mf -U)Hfi b>= _L_ J dSa J dSh A^- % ^)K{x a ,x h ,t f - tl ) . 

(4.64) 

Representing, in the usual way, the kernel by the path integral, we conclude that the 
S-matrix elements Sab are given by a Fourier transform of the path integral. Either per- 
turbative or semi-classical expansion techniques, as described earlier on in this section, 
can now be employed to obtain series expansions for these S-matrix elements. 

4.7 Exercises 

1. Using the Heisenberg picture operator equations of motion or Ehrenfest's theorem, 
show that for a (possibly time-dependent) harmonic oscillator the time-ordered 
2-point function is a Green's function, 

{-d 2 t - oj{tf) < 0\T(x(t)x{t'))\0 >= . (4.65) 

2. Using the generating functional Zo[j] (4.20), express the 4-point function 

< 0\T(x(ti)x(t 2 ))x(t 3 )x(U))\0 > (4.66) 
in terms of sums of products of Green's functions (2-point functions). 

3. Calculate the order h correction to the stationary phase approximation of the 
integral 

' / + °° dx e (WW = - % - e (*/ft)/(*c) (1 + . .) + . . .) (4.67) 

J-oo \ f"(x c ) 



y/2m% J-oo \/f"( x c 
Note that two different terms contribute to the integral at this order. 
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A Gaussian and Fresnel Integrals 

A.l Basic 1-dimensional Integrals 

The basic Gaussian integral is 

"+ 00 , _™.2/o /2tt 



/+°° 2 /O /27T 

dxe- ax / 2 = \ — q£1,«>0 . (A.l) 
-oo V a 

This result can for instance be established by the standard trick of squaring the integral 
and passing to polar coordinates. The first generalisation we will consider is the integral 

h [a, j] = [ +0 ° dx e ~ ax2 l 2 + i x jGK. (A.2) 

J — oo 

By completing the square and shifting the integration variable (using translation imparl- 
ance of the measure) one finds 



h [a, j] = r dx e " ^/«) 2 /2 + f/2a = [E e J 2 /2a (A . 3) 



Another proof of this identity is based on the trivial identity 

f+oo A 2 



f °° dx d-ax 2 /2 + jx = Q (A 4) 

Expanding this out, one finds 

= ( + °° dx (- ax + j ) e -ax 2 /2 + jx 
•j — 00 



+00 



-00 
d 

ot-Q-.h[&,j]+ jh[oi,j] ■ (A.5) 



This differential equation for I\[a,j], 

§jh[a,j] = (j/<*)h[a,j] , (A.6) 

called a Schwinger-Dyson Equation in the quantum field theory context, is evidently 
solved by 

h[a,j]= cei 2 / 2a , (A.7) 
and the normalisation (initial) condition /i[a,0] = Iq[o\ then leads to 

h[a,j]=I [ a ]ei 2 / 2a , (A.8) 
which agrees with the result (A. 3). 
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The other generalisation that will interest us is the oscillatory Fresnel integral 

/+°° • 2 In 

dx e mx I 2 . (A.9) 
-oo 

It can be obtained by noting that the basic Gaussian integral is well defined for a £ C 
with Re(a) > 0, and that it continues to be well defined for Re(a) — > provided that 
Im(a) / 0. It follows that 



J [a] = = /feW4 . (A . 10) 

A useful way to remember this result is to relate it to the integral Iq with a phase 
depending on the sign of a, 

Ma} = I [\a\]e^(^/^ . (A.11) 
By the same reasoning, for the integral 

J l [a,j]=h[-ia,ij]= [ +0 ° dxe iax2 / 2 + i i x (A.12) 

one finds 



Ma.j] - . r -l e ^J 2 /2a (A.13) 



For an alternative proof of this identity see section 4.5. 



As an application, we consider an integral of the kind we encountered in sections 1 and 
2, namely 

K = ^ f _J +0 °dpe( i / fl ^Py- e P 2 / 2m ^ . (A.14) 

By completing the square, shifting the integration variable and using the Fresnel integral 
formula, or directly upon using the above result for Ji[a,j], one finds 



K = J_ /ZTrnm (i/n)e(m/2)(y/€) J 
2vr?i V ie 



— -eh 2 e . (A.15) 

A. 2 Related Integrals 

By symmetry considerations, one has 

/+°° 2 In 

dxe -ax /2 x 2m+l = Q _ ^ A _ 16 j 
-oo 
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The integrals of even powers x 2m can be calculated by relating them to the integrals Iq 
or I\. For example, one has 

Qm /-+QO 



dxe -ax /2 x 2m = t_ 2) m° dxe -ax /2 

y ' da m J-oc 



Qva 



(A.17) 



This can be evaluated to give 

r+oo 



/+oo 2 In 

dxe -ax /2 x 2m = (_2)^V27 
-oo 



5c/ 



a 



'1/2 



J — oo 



the generating function 



= V2^(2m- l)\\a- {2m+1)/2 . (A.18) 

However, this formula does not generalise in any useful way to path integrals. An 
alternative (and very slick) way to obtain the result (A.18) is to introduce, for the 
integrals 

(A.19) 
(A.20) 

(A.21) 

(A.22) 
(A.23) 



By evaluating the sum, we see that 

"+00 



Z(j)= I °dxe-^ 2 /2 + JX = /iKj] = 

J — oo 

Hence, comparing powers of j, one deduces ^2m+i = and 



J 2 /2a 



7 l^tO^-m^)- 

Zim = \\ — (2a) j— • 

V a ml 

In view of 

(2m)! = 2 m m!(2m-l)H , 
this agrees precisely with (A.18). 



Alternatively, instead of comparing powers, one can differentiate with respect to j, so 
that the result can be written as 



/+oo 
-oo 



dx e ax2 / 2 x e = 



pft r+oo 2 in , • 

— / dxe~ ax l 2 + 3 x 



df J- 




3=0 



df 



h[a,j] 



3=0 



More generally, one has 



+ °° 2 In 

dxe~ ax l 2 F{x) 



h[a,j] 



(A.24) 



(A.25) 



3=0 
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and 



f 

J — I 



• dxe - ax y2+F( x ) 



p ( JL) 
e \ 8 j)h[a,j] 



(A.26) 

3=0 

for any function F(x). This justifies the name generating function and explains why 
the central object of interest for integrals with Gaussian weight is the integral Z(j) = 
h[a,j]. 

A. 3 Gaussian Integrals and Determinants 

The simplest 2-dimensional Gaussian integral is 

rdxTdye-™ 72-/V/2 = J^ . (A . 27) 
J-oo J-oo yap 

This generalises to 

r dx r dye-^ + & + 2 ^' 2 = , 2 " 9 , (A.28) 

J-oo J-oo \/a/3 — 'y 2 

as can be seen by completing the square. 
In terms of the matrix 

*-{*],) (A ' 29) 

the exponent of the integrand can be written as A a j ) x a x h with x a = (x, y) and the above 
identity reads 



Sx e -a^J-/2 = = ( det Ay (A 30) 

This way of writing this result generalises to d- dimensional integrals, d > 2. Let A be 
a symmetric, positive real (d x d)-matrix. Then one has 

/ + oo °°^e-^V2 = ( det _i)- 1/2 . (A .3i) 

Likewise, for the oscillatory version of this integral one has 

,d iA„,x a x b /2 ( , A Y 1 ' 2 



f 

J — i 



d d x e iA ab x-x b /2 = (^ det _A_y (A32) 



Thus the task of calculating Gaussian integrals and their various relatives and descen- 
dants has been reduced to the purely algebraic task of calculating determinants of 
matrices. 

The proof of the fundamental identities (A. 31, A. 32) is remarkably simple and will be 
left as an exercise. The strategy is to first prove them for A diagonal (this is trivial) and 
to then use the fact that any real symmetric A can be diagonalised by an orthogonal 
transformation, combined with the rotation invariance of the measure d d x. 
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B Infinite Product Identities 



A useful identity, ubiquitous in determinant and index calculations, is the following 
infinite product representation for (sinx)/x, 



00/ 2 

SIM -t-r X 



n(i-^l CB.1) 



n=l 



The validity of this formula is plausible as the right hand side has exactly the same zeros 
and identical behaviour as x — > as the left hand side. It can be proved rigorously in a 
variety of ways. The easiest (and closest in spirit to the rough argument in the previous 
sentence) is to extend (sinx)/x to a holomorphic function (sin z)/z in the complex plane 
and factorise using the Mittag-Leffler pole expansion - see e.g. section 7 of 

G. Arfken and H. Weber, Mathematical Methods for Physicists, Elsevier (2005). 

The single identity above can be used to generate an infinite number of other identities, 
in the spirit of Euler. The first non-trivial one of these results from comparing the 
quadratic term in the expansion of the left hand side, 

sinx = x - x 3 /6 + . . . = 1 _ x 2 + 
x x 6 

with the quadratic term arising from the right hand side, 

00 / 2 \ 00 2 

n=l \ / n=l 

which provides a cute proof of the famous identity 

°° 1 7T 2 

n=l 

There are a host of other occasionally useful infinite product identities (none of which, 
however, will be used in these notes). For example, closely related infinite product 
representations of other trigonomotric and hyperbolic functions are 

sinh, » ^ ^\ 

C ° SX = 11 (R6) 
cosh, = 11 [l+ (n _i / 2 2) 2 7r2 ) (B.7) 



00 

m 


(<* 


n=l 




00 

m 


(- 


n=l 




00 

m 


('• 


n=l 
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There is also the more mysterious identity 




(B.8) 



As another example, Euler's identity for the Riemann (^-function 

CM = E - s (B.9) 

n=l 

is 

oo -. 

cm = n ' (B - 10) 

n=l 1 ~ 

where p n = 2, 3, 7, 11, . . . is the sequence of prime numbers. This identity provides the 
cornerstone of the relation between number theory and complex analysis. 



C Equivalence of the VVPM and GY Formulae for the Fluctua- 
tion Determinant 



The article of Kleinert and Chervyakov cited in section 3.5 contains a nice proof of the 
classical identity (3.63) expressing the equivalence of the VVPM (3.49) and Gelfand- 
Yaglom (3.51) results for the ratio of fluctuation determinants. 

Consider the classical solution x c (t) with x c (tf^) = Xfj. This solution can equally well 
be regarded as a function of the initial position X{ and velocity Xj, 

X c — X c (t, Xi, Xi) . (CI) 

Writing this as a linear combination of two linearly independent solutions fi(t) and 
f2(t) of the oscillator equation, 

x c (t,Xi,Xi) = Xifi(t) + Xif 2 (t) (C.2) 

and imposing the conditions x c (U) = X{ and x c (tj) = Xi, one finds 

x c (U) = Xi h{U) = 1 f 2 (U) = 

x c (U) = Xi => h{U) = f 2 (U) = 1 . (C.3) 

This shows that f 2 (t) = F u (t) is the GY solution, 

F ^ t) = 9x c (t Xi,xi) (C 4) 

UXi 

At this point, one can either use directly the Hamilton- Jacobi relation 

dS Q^ = ~Pi = -mxi , (C.5) 
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to deduce (3.63). Or one can evaluate explicitly the classical action in terms of F^t) 
and fi(t) = Gu(t) (the "dual" GY solution, which plays a role analogous to the GY 
solution when one considers periodic or anti-periodic boundary conditions instead of 
zero boundary conditions). Thus one has 

x c (t) = XiG^t) + ±iF u (t) . (C.6) 

and therefore, in particular, 

x f = XiG u (tf) + XiF w (tf) 

if = XiG w {t f ) + XiF u {t f ) . (C.7) 
Since only boundary terms contribute to the classical action, it is simply given by (3.48) 

771 

s \ x c] = y Kf^c(if) -XiX c (U)] . (C.8) 

Using (C.7) to eliminate ±ij in favour of Xij, and using the fact that the Wronskian of 
F u (t) and G u (t) is t-independent, 

{F u G u -F u G u )(U) = -l = (F u G u -F u G u )(tf) , (C.9) 

one finds that the classical action, as a function of Xi and Xf, is 

S[Xc] = 2FUti) [Pu}{tf)X " + G »W*i ~ 2xiX A ■ ( C - 10 ) 



It follows that 

d 2 S[x c ] m 



(C.ll) 



dxidxf F w (tf) 
as we set out to show. 

D Poor Man's Proof of some (-function Identities 

The purpose of this appendix is to give a heuristic proof of the C-function identities 

C(0) — | (£~ii) 

C(-l) = -& ( T.n=in ) (D.l) 
C("2) = (E^=in 2 ) 



where 



oo 



C(-) = E n ~ S ■ (D-2) 



n=l 

Let us, for reasons that will become apparent below, start with the sum 



n=l 
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which we consider as the regularisation of X^Li 1 (to which it reduces as e — > 0). This 
sum is elementary, 

The derivatives of Z(e) are 



ne ~ en 



n=l 

oo 

Z"( e ) = ^n 2 e"^ , (D.5) 

n=l 

etc. Therefore, just as we formally have Z(0) = J2 n 1> we nave 

oo 

n=l 

oo 

]Tn 2 = Z"(0) , (D.6) 



n=l 



etc. Now let us expand Z(e) for small e, 



Z{e) = \-\ + ^ + O(^) (D.7) 
(in particular, there is no e 2 term). Therefore we also have 

Z"(e) = | + 0(e) . (D.8) 

Evidently, in each case the first term is singular as e — > 0. Now comes the trickery. 
Assume that, for present purposes, the regularisation (analytic continuation) of the 
Riemann ("-function amounts to nothing more and nothing less than that this term is 
absent. Then we find 

oo 2 

( X! tyreg = Z(0)reg = ~ » 
n=l Z 

oo 

(£™)re 3 = -^'(0)reg = - — 
n=l 

oo 

(E = ^"(0)re 9 = , (D.9) 

71=1 

which agrees precisely with the values of £(s) for s = 0, —1, —2 respectively. 
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As an aside: an elegant way to obtain the analytic continuation of the Riemann £- 
function is via the (assumed known) analytic continuation of the Euler Gamma-function 



f°° + 

T(s) = / dte-H 3 - 1 , (D.10) 
Jo 

which generalises the factorial, 

r(n+ 1) = n! iiGN . (D.ll) 
The first step is to establish the identity 

r(a)C(a) = / dt -— , (D.12) 
JO e 1 - 1 

which follows from changing variables t — > nt in (D.10), 

r(a) = n s r dt e'^t 3 - 1 , (D.13) 
Jo 



and 

00 roc roc fs 

rWCW = E / '//<■ 1 I 'it (! , CD.u) 



From this, the expansion of (e* — 1) 1 (as above) and the analytic continuation of the 
Gamma-function, one can then obtain the analytic continuation of the ("-function. 
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